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ABSTRACT 

We report deep Green Bank Telescope (GBT) spectroscopy in the redshifted NH3 (1,1), CS 1-0 and 
H2CO Ooo-loi lines from the z ^ 0.685 absorber towards B0218H-357. The inversion (NH3) and rotational 
(CS, H2CO) line frequencies have different dependences on the proton-electron mass ratio /i, implying that a 
comparison between the line redshifts is sensitive to changes in /i. A joint 3-component fit to the NH3, CS, 
and H2CO lines yields [A/i//i] = (—3.5 ± 1.2) x 10^^, from z ^ 0.685 to today, where the error includes 
systematic effects from comparing lines from different species and possible frequency-dependent source mor- 
phology. Two additional sources of systematic error remain, due to time variability in the source morphology 
and velocity offsets between nitrogen-bearing and carbon-bearing species. We find no statistically-significant 
(> 3a) evidence for changes in /i, and obtain the stringent 3cr constraint, [A/i//i] < 3.6 x 10~^, over 6.2 Gyrs; 
this is the best present limit on temporal changes in /i from any technique, and for any lookback time, by a 
factor > 5. 

Subject headings: atomic processes — galaxies: high-redshift — quasars: absorption lines 



1. INTRODUCTION 

A generic prediction of theories that attempt to unify the 
standard model of particle physics with general relativity is 
that particle masses and low-energy coupling constants vary 
with space and time. A detection of such changes would im- 
ply new physics beyond the standard model. Tests of changes 
in fundamental "constants" such as the fine structure constant 
a, the proton-electron mass rat io u or the p roton g-factor gp 
are hence of much interest (e.g. IUzanll2003l) . 

Redshifted atomic and molecular spectral lines provide an 
interesting avenue to probe the possibility of fundamental 
constant evolution over cosmological timescales (Savedoff 
Il956h . Comparisons between the redshifts of lines whose 
rest frequencies have different dependences on a con- 
stant like a are sensitive to changes in the constant [see 
iKanekaij (l2008h for a recent review]. There are a num- 
ber of such technique s, using cornparis o ns between differ- 
ent transitions (e.g. Thompson 1975; 'Wolfe et al. 1976 



Dzubaetal. 1999; Darling 2003; Chengalur & Kanekar 2003, 
Kaneka r & Chengalur 2004; Flambaum & Kozlov 200^. In- 



deed, one of these techniques, the many- multiplet method , 
has yielded evidence for changes in a: Mur phy et all (120041) 
obtained [Aa/a] = (-5.4 ± 1.2) x 10"''^ from 143 ab- 
sorbers with (z) — 1.75 (see also Murphy et al. 2003). 
This result has as yet been neither confirmed nor ruled out 
(e.g. 'Srianand et aLll2007l iMolaro et al.ll200l iMurphy et al.] 
|2008b; KanekEjetallllQlOb), but evidence has now been 
found for additional systematic effects in the data, due to er- 
rors in the wavelength calibration (Griestetal. 2010). Re- 
cently, an indepe ndent te chnique, using radio "conjugate" 
satelUte OH lines (IKanekar et al.,.200 4; Kanekar et al. 2005.) . 
found weak evidence (at ~ 99.1% confide nce level) for 
changes in a combination of a, fi and gp: IKanekar et al.l 
(|20 10a) obtained [AG/G] = (-1.18 ±0.46) x IQ-^ between 
z ^ 0.247 and the present epoch, where G = gp[^a'^]^'^^ . If 
changes in a are assumed to dominate over those in /i and gp, 
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one obtains [Aa /a] = (—3.1 ± 1.2) x 10~^, consistent with 
the result of Mu rphy et al.l (120041) , albeit at a lower redshift. 

For three decades, ro-vibrational molecular hydrogen (H2) 
lines provided the sole method to directly probe changes in 
the proton-electron mass ratio (iThompsonI 1 1 9751) . Unfortu- 
nately, f ew redshifted H2 ab sorbers have been detected so 
far (e.g. iLedoux et al.ll2003b , and, of these, only four have 
been found suitable to study changes in /i (e.g. Ivanchik et aLl 
"2005; Reinhold et al. 2006; King et al. 2008; Thompson et al] 
2009; Malec et al.l 120101: LWendt & Molaro 2011). These re- 
sults have also been controversial: Reinhold et al.^ (.200 6) ob- 
tained [A^/^] = (+2.0 ± 0.6) x 10~^ from two absorbers 
at z ^ 2.6 — 3.0, i.e. evidence for a larger value of ^ 
at high redshift (see also Ivanchik et al. 2005). However, 
independen t re-analy ses of these data by King et al. (20081) 
and Thompson et al.l (i2009), using improved wavelength cal- 
ibration techniqu es, found no evidence for changes in /i. 
iKing et al.l(l2008l) obtain [A/i//z] = (+2.6±3.0) x lO'^ from 
th ree absorbers at z ~ 2.6 — 3.0 (including t he two systems 
oflReinhold et al.l200 6^); note, however, that Wendt & Molaroj 
(Ebii) argue that systematic errors have been under-estimated 
in the latter analysis. 

A new l aboratory technique to prob e changes in /i was pro- 
posed by Ivan Veldhoven et"an (l2004b , usmg mversion tran- 
sitions of deuterated ammonia (ND3). 'Flambaum & Kozlov| 
^)07) adapted this technique to astrophysical circumstances, 
using ammonia (NH3) inversion lines and rotational lines. 
Rotational and inversion line frequencies have different de- 
pendences on /i, and a comparison between the redshifts of 
NH3 inversion lines and rotational (e.g. CO, HCO+, CS, etc) 
lines from a single cosmologically-distant absorber is thus 
sensitive to changes in /i. Only two redshifted NH3 absorbers 
are currently known, at z 0.685 towards B02 18-1-357 
(Eenkel et al. 2005) and z - 0.886 towards B1830-210 
dHenkel et al.n2008i) . In both cases, the NH3 detection spec- 
tra have been used to obtain initial constra ints on changes in ^ 
dMurphy et al."2008a;'Henkel etalj2009i). Finally, the ammo- 
nia t echnique has also been used to search for spatial changes 
in fj,: iLevshakov et al.l(l2010l) obtained the (conservative) con- 
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straint [A/i//i] < 3 x 10^^ from studies of multiple molecular 
clouds in the Galaxy. 

The z ^ 0.685 absorber towards B0218+357 has a far 
smaller velocity sp read than the z ^ 0.886 syst em to- 
wards B1830-210 ( Wik Und & Combes|[T99 5'. 1996b), mak- 
ing it a better candidate for accurate redshift measure- 
ments. The carbon monosulfide (CS) 1-0 line [rest frequency: 
48.990957(2) GHz; iKim & Yamamotoll2003ll is a good rota- 
tional transition for the comparison with the NH3 inversion 
lines because it is known to be unsaturated in the z ^ 0.685 
absorber (Combes et al. 1997) and it lies at a frequency rela- 
tively close to that of the NH3 (1,1) lines. The formaldehyde 
(H2CO) Onn-lni line [rest fr equency: 72.837948(10) GHz; 
ICornet & Winnewisserlll980ll was chosen as a second rota- 
tional transition, to obtain an estimate of systematic effects 
in the above comparison; th is too is known to b e unsaturated 
in the z ^ 0.685 absorber ("Jethava et al."2007'). This Letter 
reports Green Bank Telescope (GBT) spectroscopy in the red- 
shifted NH3 (1,1), CS 1-0 and H2CO Ooo-loi lines, that yield 
the best present constraints on changes in the proton-electron 
mass ratio. 

2. OBSERVATIONS, DATA ANALYSIS AND RESULTS 

The GBT observations of the NH3 (1,1) and CS 1-0 lines 
at z ^ 0.685 towards B0218H-357 were carried out in January 
and August 2008 (proposals AGBT07C-016 and AGBT07C- 
054), using the Ku- and Ka-band receivers, respectively. 
A spectrum in the redshifted H2CO Oqo-Ioi line was ob- 
tained in July 2008 with the GBT Q-band receiver, as part 
of an absorption survey of the z ^ 0.685 galaxy (pro- 
posal AGBT07C-054). Bandpass calibration was carried out 
via beam-switching for the NH3 observations, and with sub- 
reflector nodding for the higher-frequency CS and H2CO ob- 
servations. Two polarizations were used for the observations 
of the NH3 (Ku-band) and H2CO (Q-band) lines, and a sin- 
gle Ka-band polarization for the CS line. The total observ- 
ing times were 6 hours (NH3), 1.5 hours (CS) and 3 hours 
(H2CO), with the NH3 observations split into two observing 
sessions, in January and August 2008. Standard procedures 
in the packages DISH (for the NH3 Hnes) and GBTIDL (for 
the CS and H2CO lines) were used to analyse the data. The 
final NH3, CS, and H2CO spectra are shown in Figs. [IlA], 
[B] and [C]; these have velocity resolutions of ^ 0.26, 0.12 
and 2.8 km/s, respectively (after Hanning-smoothing and re- 
sampling), and root-mean-square noise values of ~ 0.00082, 
0.0097 and 0.0026, in optical depth units, per independent 
velocity channel. 

A GBT Ku-band search was also carried out for the CCS 1- 
transition at z ~ 0.685; no absorption was detected at the 
redshifted CCS line frequency, with a 3a optical depth limit 
of ^ 0.0016 per 4.2 km/s channel. The K- and Ku-band re- 
ceivers were further used to search for redshifted NH3 (1,1) 
absorption from two other redshifted molecular absor bers, at 
z - 0.247 towards B1413H-135 (Wiklind & Combes '199^, 
and z - 0.674 towards B1504H-377 (Wiklind & Combes 
Il996ah . with no detected absorption in either case. The 
3(T optical depth limits were 0.00068 per 4 km/s channel 
(B1413H-135) and 0.00061 per 4 km/s channel (B1504H-377). 

The test for changes in the proton-electron mass ratio fi was 
carried out through a simultaneous fit to the NH3 (1,1), CS 1- 
and H2CO Ooo-loi spectra, using the package VPFIlQ. The 
velocity structure in all optically-thin lines was assumed to 

^ http://www.ast.cam.ac.uk/ffc/vpfit.html 



be the same, in both the number of absorbing "clouds" and 
their redshifts. Turbulence was assumed to be the dominant 
contributor to the line widths, which were hence tied together 
in the fit (Murphy et al. 2008a). The NH3 hyperfine struc- 
ture was included by assuming local thermodynamic equi- 
lib rium (LTE), with the hyperfine ratios taken from Table S 1 
of (iMurphv et al.ll20(j8ah . Finally, the fit also included a sin- 
gle velocity offset between the spectral components in the in- 
version and rotational transitions, to account for a possible 
change in the proton-electron mass ratio. 

A simultaneous multi-component Voigt profile fit was car- 
ried out to the three spectra of Fig.[T|with the above assump- 
tions, aiming to minimize xt varying the fit parameters. 
A 3-component fit was found to yield xt — 105 and no ev- 
idence for structure in the residuals after subtracting out the 
fits to each spectrum (see the lower panels of Fig. Specif- 
ically, a Kolmogorov-Smirnov rank-1 test found the fit resid- 
uals for all lines to be consistent with a normal distribution 
within ^ la confidence. One- and two-component fits were 
found to yield both a larger and clear non-Gaussian struc- 
ture in the residual spectra. While increasing the number of 
spectral components to four did yield a marginally-lower 
(~ 1.02), the additional spectral component was only weakly 
detected in the H2CO spectrum (at < 3(t significance) and 
was not visible in the NH3 or CS transitions. A 3-component 
model thus appears to provide a good fit to the data; however, 
for completeness, results from the 4-component fit will also 
be mentioned below. 

For the coupled 3-component model, the number of fit pa- 
rameters is 16, three redshifts, three line widths, nine optical 
depths in the NH3, CS, and H2CO lines, and the velocity off- 
set between the inversion and rotational transitions. The pa- 
rameters of the best 3-component fit are listed in Table 1; the 
best-fit velocity offset is AV^ = (-0.36±0. 10) km/s, with the 
NH3 lines blueshifted relative to the rotational lines. A sim- 
ilar result is obtained from the best-fit 4-component model, 
AV^ ~ (—0.39 ±0.11) km/s. Similar results are also obtained 
from carrying out 3-component fits to the NH3 and CS lines 
(i.e. without including the H2CO line), and to the NH3 and 
H2CO lines (i.e. without the CS Une). 

To obtain an estimate of systematic effects, a similar 3- 
component fit was also carried out to the CS and H2CO ro- 
tational lines alone (i.e. without including the NH3 inver- 
sion lines). This yielded a velocity offset of 6V = (0.029 ± 
0.068) km/s, between the CS and H2CO lines. The la error 
obtained here, when comparing two rotational lines, will be 
used as an estimate of the systematic error due to local veloc- 
ity offsets between different species in the absorber 

The fractional change in the proton-electron mass 
ratio /i is related to the measured velocity offset 
between inversion and rotational line s AT^ by the 
expression (Flambaum & Kozlov l2007h [A^/^] ~ 
0.289 [(z,™ -z™t)/(l + z)] K. 0.289 [A V'/c], where 
Zinv and Zrot are, respectively, the inversion and rotational 
redshifts, and z is their average. Our final velocity offset is 
AF = [-0.36 ± 0.10(stat.) ± 0.068(syst.)] km/s; this yields 
[A/x/^x] = [-3.47 ± QM{stat.) ± {)M{syst}}\ x 10"^. 
Adding the statistical and systematic errors in quadrature, we 
obtain [A^/^] = (-3.5 ± 1.2) x lO"'^, between z - 0.685 
and the present epoch. 

3. DISCUSSION 

An excellent summary of the systematic effects inherent 
in the comparison between rotational and inversion lines in 
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Fig. 1 .— GBT spectra in the redshifted [A] NH3 (1,1), [B] CS 1-0 and [C] H2CO Ooo-loi transitions from the z ~ 0.685 absorber towards B0218+357. The 
upper panels show normalized line depth in each transition plotted against velocity, in km/s, relative to a heliocentric redshift of z = 0.684676; the 3-component 
fit is overlaid on each spectrum. The lower panels show the residuals from the spectra after subtracting out the joint 3-component fit; these are consistent with 



the z 0.685 absorber is given by iMurphv et alj (l2008al) . 
These include the following: (1) the assumption that dif- 
ferent transitions have the same number of spectral compo- 
nents, and at the same redshifts, (2) the background source 
morphology is frequency-dependent, implying that the NH3 
and rotational transitions might arise along slightly different 
sightlines, (3) the background source flux density (and hence 
morphology) varies with time, so observations at different 
epochs might probe different sightlines, (4) the assumption of 
LTE for the NH3 hyperfine structure, (5) saturation effects if 
highly-saturated lines [e.g. the HCO"*" and HCN lines used by 
iMurphy et alJ (|2P08a)] are used, and (6) local velocity offsets 
be tween the species giv ing rise to the transitions. As noted 
bv lMurphy et alJ (l2008al) . all these effects could yield signif- 
icant contributions to their systematic eiTors. For example, 
the inversion and rotational transitions used by lMurphy et alJ 



(l2008a l) were at very different frequencies, ^ 14 GHz and 
^ 106 GHz; changes in the background source morphology 
could thus result in different sightlines in the different tran- 
sitions. The NH3 and HCO+/HCN observations were sepa- 
rated by a few years, implying that temporal variability in the 
source morphology could also be an issue, especially given 
that saturated transitions were used in the analysis. Finally, 
the NH3 spectra were of too low signal-to-noise ratio (S/N) 
to detect the NH3 hyperfine structure and to test for non-LTE 
effects. 

Most of the above issues have been directly addressed in 
the present work. All three transitions used in the analy- 
sis are unsaturated and have been observed at high S/N, and 
with the same telescope. The two outlying NH3 hyperfine 
components (at ±19.5 km/s relative to the strongest compo- 
nent) in the NH3 (1,1) line have also been clearly detected 
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Component 


Velocity offset 


Redshift 


FWHM 


Line depth X 100 






km/s 




km/s 


NH3 


CS 


Ha CO 


1 




0.6846935(2) 


3.99 ± 0.13 


0.496 ±0.060 


13.83 ±0.78 


15.61 ±0.84 


2 


(-0.36 ±0.10) 


0.6846858(10) 


10.60 ±0.45 


0.734 ±0.026 


4.39 ±0.33 


3.00 ±0.39 


3 


0.6846214(7) 


3.53 ±0.33 


0.128 ±0.032 


2.84 ±0.28 


5.17 ±0.37 



TABLE 1 

Parameters of the best 3-component Voigt profile fit to the NH3 ( 1 , 1 ), CS 1-0 and H2CO Goo- 1 01 spectra. The second column 

CONTAINS THE BEST-FIT VELOCITY OFFSET BETWEEN THE ROTATIONAL AND INVERSION LINES; NOTE THAT A SINGLE VELOCITY OFFSET WAS 
ASSUMED FOR ALL SPECTRAL COMPONENTS . THE REDSHIFTS OF COLUMN (3) ARE IN THE HELIOCENTRIC FRAME, WHILE THE LINE DEPTHS (XlOO), 
LISTED IN THE LAST THREE COLUMNS, HAVE BEEN NORMALIZED BY THE SOURCE CONTINUUM AT EACH OBSERVING FREQUENCY. 



(see Fig. [HA]). As noted by lMurphy eTaP (l2008ah . non-LTE 
effects resulting in hyperfine "anomalies" should cause the 
satellite hyperfine components (with AFi ^ 0, where F\ 
is the quadrupole quantum number) to have different optical 
depths. No evidence for such non-LTE conditions is apparent 
in the ±19.5 km/s satellite hyperfine components in Fig.[TtA]; 
the optical depths are found to agree within the noise. 

We note, finally, that the systematic error of 7.6 x 10^^ 
in the result of Murphv et al. (2008a) contains two contribu- 
tions: (1) 7 X 10^^ because spectral components detected in 
the HCO+ and HCN spectra may not be detected in the NH3 
spectra, due either to the large frequency difference between 
the inversion and rotation lines or to the fact that the HCO+ 
and HCN lines are optically thick, and (2) 3 x 10~^ due to 
the possibility of non-LTE effects in the NH3 spectra. Our 
choice of rotational lines and the higher sensitivity of our NH3 
spectra implies that neither of these are significant sources of 
s ystematic error in t he present analysis. 

iLiszt et al] (I2006h find that NH3 column densities correlate 
best with CS and H2CO column densities in Galactic diffuse 
clouds; CS and H2CO are thus likely to be the best rotational 
transitions for the inversion/rotation comparison, as the corre- 
lation in column densities suggests that the three species are 
likely to arise in the same part of a gas cloud. Galactic absorp- 
tion has also been detected in all three species to wards two 
quasars, B0355-H508 and B0415-H379 (ILiszt & L ucas 1995; 
iLucas & Liszt"1 998l: ILiszt et all 120061) ; both sightHnes show 
two components in all specie s, with line velociti es agreeing 
within 0.2 km/ s [see Table 1 of lLiszt & LucasI 1 995t Table 6 of 
ILucas & Lisz3[l998l and Table Al of lLiszt et al.l l20061; local 
velocity offsets between the species thus appear to be small. 
Further, the rotational and inversion line frequencies are much 
closer here than in the comparison of Murphv et al. (2008a), 
with the redshifted NH3, CS, and H2CO lines at - 14 GHz, 
~ 29 GHz and ^ 43 GHz, respectively. Changes in source 
morphology with frequenc y are thus a le ss important issue 
here than in the analysis of 'Murp hv et al.l (l2008a) . A direct 
test comes from the comparison between the two rotational 
lines, whose frequencies differ by a factor of ^ 1.5, compa- 
rable to the ratio of the CS and NH3 frequencies 2). The 
CS-H2CO comparison also provides an estimate of the sys- 
tematic error due to local velocity offsets between different 
species in the clouds. This test found the CS and H2CO line 
redshifts to agree within the noise; the Ict error in this com- 
parison has been used to quantify the systematic error due to 
both local velocity offsets and differing background source 
morphology at the different line frequencies. Finally, similar 
velocity offsets (consistent within the errors) were obtained in 
the independent comparisons between the NH3 and CS lines, 
and the NH3 and H2CO Hnes. 

A source of systematic effects that could not be directly 



addressed here is time variability in the background source 
morphology, which might yield different sig htUnes through 
the absorbing clouds at different epochs [Murphy et al. 2008at 
see Muller & Guehn (2008) for the z ~ 0.886 lens towards 
B1830-210]. Unfortunately, the redshifted NH3, CS and 
H2CO transitions towards B0218H-357 require different GBT 
receivers and cannot be observed simultaneously. The present 
NH3, CS and H2CO observations were carried out between 
January and August 2008; the possibility that the weak offset 
between the inversion and rotational lines arises due to small 
changes in the sightline thus cannot be ruled out. Further, the 
agreement between CS and H2CO velocities does not rule out 
the possibility that nitrogen-bearing species like NH3 arise at 
different velocities than carbon-bearing species like CS and 
H2CO. This can only be tested by using rotational transitions 
of other nitrogen-bearing species (e.g. HC3N, CH3CN, etc). 

While there are two possible sources of systematic error 
that we are as yet unable to quantify, the present data show 
no statistically-significant (> 3(j) evidence for changes in the 
proton-electron mass ratio. Our 3cr upper limit on changes in 
/i is [A/i//i] < 3.6 X 10^'', between z ^ 0.685 and z = 
(a lookback time of 6.2 Gyrs). For comparison. iMurphv et all 
(i2008a ) obtained [A/i//i] < 1.8 x lO ''^ (2g) in the z ~ .685 
absorber towards B0218H-357, while iHenkel et all (120091) ob- 
tained [A^/^i] < 1.4 X 10-^ (3cr) in the z ~ 0.886 ab- 
sorber towards B1830— 210, with both results b ased on inver- 
sion/rotation comparisons. Note, however, that iHenkel et al.l 
(2009) used single-Gaussian fits for the NH3 and HC3N lines 
towards B1830— 210, although it is clear from the profiles of 
other unsaturated lines (e.g. CS) that at least three absorb- 
ing "clouds" are present along the sightline. Further, the GBT 
spectra in the strongest [(1,1), (2,2) and (3,3)] NH3 lines to- 
wards B1830— 210 were severely affected by radio frequency 
interference (see Fig. 1 of iHenkel et al. 2008). The effect 
of these issues on the results of iHenkel et al.i (i2009i) is un- 
clear, but the error budget is likely to increase. At higher 
redshifts, the best published constraint on changes in /i is 
that of King et al. (2008) : [A/x/fx] < 6.0 x 10"^ using H2 
lines (but see Wendt & Molarol [201 ll) . The present result, 
[A/x//x] < 3.6 X lO"'^ (3cr) at z ~ 0.685, is thus the most 
sensitive constraint on temporal changes in /i at any redshift, 
by a factor > 5. 

To compare this result with those from laboratory stud- 
ies, it is necessary to assume a model for the variation of 
[I with time. For linear variation, the present result yields 
[ij [I < 5.6 X 10^^^ per year, more than an order of magni- 
tude better than the be st (model-dependent) la boratory con- 
straint on changes in (iRosenband et al.ir2008l) . 

In conclusion, we have used the GBT to obtain high 
S/N spectra in the NH3 (1,1), CS 1-0 and H2CO Ooo- 
loi transitions from the z ^ 0.685 absorber towards 
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B02 18+357. A comparison between the redshifts of the in- 
version and rotational lines yields a velocity offset of AV ~ 
[-0M±0.10{stat.)±0m8{syst.)] km/s, with the NH3 lines 
blueshifted relative to the rotational ones. Two sources of sys- 
tematic error remain to be quantified, arising from ( 1 ) velocity 
offsets between nitrogen-bearing and carbon-bearing species, 
and (2) time variability in the background source morphology, 
which might yield different sightlines at different epochs. We 
find no statistically-significant (> 3(t) evidence for changes in 
/i, with the 3(7 constraint [A/u//u] < 3.6 x 10~^ over a look- 
back time of 6.2 Gyrs. This is the strongest present constraint 
on temporal changes in the proton-electron mass ratio. 
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